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ABSTRACT 


'Cauchy  problems  for  equations  modelling  non-Newtonian  fluids  are 
discussed  and  recent  existence  theorems  for  classical  solutions,  based  on 

semigroup  methods,  are  presented.  Such  existence  results  depend  in  a  crucial 

y  S' 

manner  on  the  symbol  of  the  leading  differential  operator.  Both  ^parabolic* 


anc 


»d^~hyperbolic"  cases  are  discussed.  In  general,  however,  the  leading 
differential  operator  may  be  of  non-integral  order,  arising  from  convolution 
with  a  singular  kernel.  This  has  interesting  implications  concerning  the 
propagation  of  singularities.  In  particular,  there  are  cases  where 

£  -smoothing  coexists  with  finite  wave  speeds. 
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SIGNIFICANCE  AND  EXPLANATION 


This  paper  reviews  recent  existence  results  for  initial  value  problsms 
for  equations  modeling  viscoelastic  fluids.  Motivated  by  the  linear  case*  the 
equations  describing  thesa  fluids  ars  classified  as  parabolic,  hyperbolic  or 
intermediate. 

The  behaviour  of  these  different  types  in  linear  wave  propagation  is 
discussed.  In  the  hyperbolic  case,  propagation  of  shocks  occurs,  in  the 
parabolic  case  the  wave  speed  is  infinite  and  singularities  are  smoothed 
out.  The  intermediate  cases  allow  intaresting  situations  where  the  wave  speed 
is  finite,  but  smoothing  of  singularities  does  occur. 

Mathematical  techniques  applicable  to  the  solution  of  initial  value  are 
discussed  for  both  the  parabolic  and  the  hyperbolic  case.  In  both  cases,  we 
discuss  three-dimensional  problems  including  treatment  of  the  incompressi¬ 
bility  condition.  The  intermediate  case  remain  an  interesting  open  problem. 
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INITIAL  VALUE  PROBLEMS  FOR  VISCOELASTIC  LIQUIDS 


M.  Renardy 

1.  Introduction 

Viscoelastic  liquids  are  characterised  by  constitutive  laws  allowing  the 
stress  to  depend  on  the  history  of  the  deformation  [3],  [9],  [17],  [18].  This 
dependence  is  usually  assumed  to  be  local,  i.e.  the  stress  at  the  location  of 
a  given  fluid  particle  depends  only  on  the  history  of  the  deformation  gradient 
at  this  same  particle.  Following  Noll  [17],  such  fluids  are  called  "simple". 
The  constitutive  law  is  further  restricted  by  the  assumption  that  rigid  body 
motions  do  not  contribute  to  the  stress  (called  the  principle  of  frame- 
indifference)  and  by  material  symmetries.  In  this  paper,  I  shall  deal  only 
with  isotropic,  incompressible  materials. 

For  a  mathematical  existence  theory,  more  needs  to  be  known  about  the 
nature  of  the  constitutive  law.  However,  the  precise  form  is  not  known  for 
any  particular  material,  and  we  must  resort  to  models.  Rheological  models 
have  been  motivated  by  one  or  a  combination  of  the  following  considerations: 

1)  Modification  of  linear  theories  to  make  them  comply  with  frame- 
indifference. 

2)  Formal  analogy  with  finite  elasticity. 

3)  Kinetic  theories  of  chain  molecules. 

From  a  mathematical  point  of  view,  the  type  of  the  equation,  or,  in  other 
words,  the  nature  of  the  leading  differential  operator  is  of  particular 
interest.  The  rheological  models  can  essentially  be  classified  into  three 
categories : 


Sponsored  by  the  United  States  Army  under  Contract  No.  DAAG29-80-C-0041.  This 
material  is  based  upon  work  supported  by  the  National  Science  Foundation  under 
Grant  No.  MCS-82 10950. 


hyperbolic  cut<  The  equations  are  transforned  in  such  a  way  that  they  fit 
into  the  framework  of  known  theorems  due  to  Sobolevskii  [28]  and  Xato  [11-13], 
respectively.  In  both  cases  we  deal  with  three-dimensional  problems, 
including  treatment  of  the  incompressibility  condition. 


2.  A  problem  in  linear  viscoelasticity 

If  the  stress  is  a  linear  functional  of  the  strain  history,  we  can 
formally  write  the  constitutive  law  in  the  form 

t ( t )  -  a(t-s)Y(t,s)ds  .  (2.1) 

Here  t  is  the  stress  tensor  and  y  is  the  gradient  of  the  relative  dis¬ 
placement  (from  the  position  at  time  s  to  the  position  at  time  t).  In 
particular,  contributions  resulting  from  different  times  s  superpose  in  an 
additive  fashion.  Boltzmann  [3]  suggested  that  the  following  restrictions 
should  hold: 

1)  If  the  relative  strain  is  always  positive,  so  is  the  stress. 

2)  The  strain  from  a  more  remote  time  always  has  a  lesser  influ¬ 
ence  than  that  from  a  more  recent  time. 

*  A 

This  means  that  a  >  0,  a*  <  0  in  the  sense  of  distributions.  In  this  case, 
it  can  be  shown  [25]  that 

a(T)  -  -u5'(t)  +  a(T)  ,  (2.2) 

where  p  >  0,  and  a  is  a  non-negative,  non-decreasing  function  of  t  >  0. 
This  suggests  that  we  classify  constitutive  laws  according  to  the  degree 

A 

of  the  singularity  of  a  at  t  -  0.  The  strongest  possible  singularity  is 
the  4* -distribution.  If  this  term  is  absent,  it  is  important  whether  a  is 


finite  or  infinite  at  0 


We  shall  consider  the  following  particular  problem:  A  fluid  filling  the 
half-space  above  an  infinite  plate  is  at  rest  for  t  <  0.  At  t  =  0,  the 
plat  is  suddenly  set  into  uniform  motion.  We  want  to  determine  the  motion 
for  t  >  0. 

For  a  linear  viscoelastic  fluid,  this  leads  to  the  equation 
ufct(x,t)  =  puxxt<x,t)  +  a(t-s)  (uxx(x,t) 

-  u  (x,s))ds,  x  >  0,  t  >  0  , 

*  (2.3) 

u(x,t)  =0  ,  t  <  0  , 
u(0,t)  =1  ,  t  >  0  . 

In  addition  to  Boltzmann's  restrictions,  we  assume  that  a  decays  fast  enough 
at  •  so  that  the  integral  converges  (in  all  rheological  models  I  know  of  it 
decays  exponentially). 

If  w  ^  0,  the  term  uu  is  the  highest  order  term  on  the  right  hand 

xxt 

side,  and  we  would  classify  the  equation  as  "parabolic".  If  u  *  0  and  a 
is  continuous,  then  the  highest  order  tem  is  uxx  •  a(T)dx,  since  the 
convolution  is  a  differential  operator  of  order  -1.  We  would  therefore  call 
the  equation  "hyperbolic".  Intermediate  types  arise  from  singular  integral 
kernels,  which  act  like  fractional  derivatives. 

It  is  interesting  to  consider  how  these  different  types  propagate 
singularities.  Coleman  and  Gurtin  [4]  have  shown  (for  u  *  0)  that,  if  there 
is  a  propagating  shock  front,  its  wave  speed  is  /a  and  its  amplitude  is 
e”a(°)t/2A.  Here  A  =  /"  a(T)dx. 

It  was  shown  in  (16],  [24]  that  indeed  the  singularity  propagates  as  a 
shock  if  u  =  0  and  a  has  1  +  e  derivatives  in  L 1 .  For  li  t  0,  the 
shock  is  smoothed  out  and  analytic  solutions  are  obtained. 

A  class  of  kernels  which  act  like  fractional  derivatives  is  given  by 
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2.  A  problem  in  linear  viscoelasticity 

If  the  stress  is  a  linear  functional  of  the  strain  history,  we  can 
formally  write  the  constitutive  law  in  the  form 

T(t)  ■  a(t-s)y(t,s)ds  .  (2.1) 

Here  T  is  the  stress  tensor  and  y  is  the  gradient  of  the  relative  dis- 
placement  (from  the  position  at  time  s  to  the  position  at  time  t).  In 
particular,  contributions  resulting  from  different  times  s  superpose  in  an 
additive  fashion.  Boltzmann  [3]  suggested  that  the  following  restrictions 
should  holdt 

1)  Xf  the  relative  strain  is  always  positive,  so  is  the  stress. 

2)  The  strain  from  a  more  remote  time  always  has  a  lesser  influ¬ 
ence  than  that  from  a  more  recent  time. 

*  A 

This  means  that  a  >  0,  a'  <  0  in  the  sense  of  distributions.  In  this  case, 
it  can  be  shown  [25]  that 

a(T )  -  -gS*(T)  ♦  a(T)  ,  (2.2) 

where  |i  >  0,  and  a  is  a  non-negative,  non-decreasing  function  of  t  >  0. 
This  suggests  that  we  classify  constitutive  laws  according  to  the  degree 
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of  the  singularity  of  a  at  t  ■  0.  The  strongest  possible  singularity  is 
the  6* -distribution.  Xf  this  term  is  absent,  it  is  important  whether  a  is 
finite  or  infinite  at  0. 


We  shall  consider  the  following  particular  problems  A  fluid  filling  the 
half-space  above  an  infinite  plate  is  at  rest  for  t  <  0.  At  t  *  0,  the 
plat  is  suddenly  set  into  uniform  motion*  We  want  to  determine  the  motion 
for  t  >  0. 

For  a  linear  viscoelastic  fluid,  this  leads  to  the  equation 
utfc<x,t)  -  uuxxt<x,t)  +  a(t-s)(uxx<x,t) 

-  u  (x,s))ds,  x  >  0,  t  >  0  , 

xx  (2.3) 

u(x,t)  =0  ,  t  <  0  , 
u(0,t)  =1  ,  t  >  0  . 

In  addition  to  Boltzmann's  restrictions,  we  assume  that  a  decays  fast  enough 
at  “  so  that  the  integral  converges  (in  all  rheological  models  I  know  of  it 
decays  exponentially). 

If  p  f  0,  the  term  pu  is  the  highest  order  term  on  the  right  hand 
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side,  and  we  would  classify  the  equation  as  "parabolic".  If  u  *  0  and  a 
is  continuous,  then  the  highest  order  term  is  uxx  •  a(r)dx,  since  the 
convolution  is  a  differential  operator  of  order  -1.  We  would  therefore  call 
the  equation  "hyperbolic".  Intermediate  types  arise  from  singular  integral 
kernels,  which  act  like  fractional  derivatives. 

It  is  interesting  to  consider  how  these  different  types  propagate 
singularities.  Coleman  and  Gurtin  [4]  have  shown  (for  M  =  0)  that,  if  there 
is  a  propagating  shock  front,  its  wave  speed  is  /a  and  its  amplitude  is 
e  Here  A  “  J0  a(T)dr. 

It  was  shown  in  [16],  [24]  that  indeed  the  singularity  propagates  as  a 
shock  if  v  =  0  and  a  has  1  +  e  derivatives  in  L1.  For  p  t  0,  the 
shock  is  smoothed  out  and  analytic  solutions  are  obtained. 

A  class  of  kernels  which  act  like  fractional  derivatives  is  given  by 


Kernels  of  this  nature  are  suggested  by  some  molecular  theories  [7],  [26] 


[31].  As  T  ♦  0,  a  behaves  like  t  1^ot.  It  is  thus  integrable  for  a  >  1 
and  not  integrable  for  a  <  1.  Thus  the  result  of  Coleman  and  Gurtin  suggests 
an  infinite  wave  speed  for  a  <  1,  while  for  a  >  1  you  expect  a  wave  with 
finite  speed  but  with  no  amplitude.  This  is  in  fact  the  case:  If  a  >  1, 
the  solution  u  is  zero  for  x  >  /a  t  and  is  analytic  and  not  identically 
zero  for  x  <  /a  t.  Across  the  line  x  **  / A  t ,  however,  there  is  no 
singularity  and  the  solution  is  C*  [24],  If  a  <  1,  the  solution  is 
analytic  except  at  t  “  0. 

From  the  point  of  view  of  classification,  we  see  that  the  hyperbolic  term 
is  still  the  leading  order  for  a  >  1,  however,  the  order  of  the  correction 
is  lower  only  by  a  fraction,  not  by  one.  If  a  <  1,  the  right  hand  side  of 
(2.3)  is  of  fractional  order. 


3.  An  example  of  a  parabolic  model 

The  simplest  class  of  "parabolic"  models  arises  when  a  term  of  lower 
differential  order  is  added  to  a  Newtonian  term.  This  is  a  common  way  of 
modeling  dilute  polymer  solutions.  Here  you  assume  that  the  stress  consists 
of  a  Newtonian  part  arising  from  the  solvent,  and  some  additional  term  coming 
from  the  dissolved  polymer. 

Before  formulating  the  equations,  we  have  to  introduce  some  notation. 

Since  the  "simple  fluid"  is  essentially  a  Lagrangian  concept,  it  is  natural  to 

12  3 

use  Lagrangian  coordinates.  We  denote  those  by  ?  *  (C  ,C  )♦  £  varies 

over  a  bounded  domain  fl  with  a  C^-boundary.  By  we  denote  the 

position  of  the  particle  £  at  time  t.  The  deformation  gradient  F  has 
i  3y* 

components  F,  *  ~=~r»  It  has  to  satisfy  the  incompressibility  condition 
3  H3 
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I 


det  F  -  1  .  (3.1) 

T 

The  Cauchy  strain  tensor  is  defined  by  Y  *  F  F.  He  write  the  constitutive 

law  in  terms  of  the  upper  convected  stress  w,  which  is  related  to  the  Cauchy 

T 

stress  T  by  T  *  FirF  .  we  assume  that  the  constitutive  law  has  the  form 

»  =  -PY-1  "  n  < Y_1 )  +  F(y)  .  (3.2) 

A 

F  is  a  tensor-valued  functional  of  Y»  the  history  of  y  * 

A 

Y(£,t)(s)  =  Y(£,t+s),  s  e  (—,01  .  (3.3) 

If  F  «  0,  [3.2]  describes  a  Newtonian  fluid. 

The  equation  of  motion  reads 

py1  »  JL  (iL-  *rs)  +  ,  (3.4) 

3?  H 

and  we  impose  Dirichlet  conditions  on  the  boundary 

yi(£,t)  *  ♦i(i,t),  i  e  an  .  (3.5) 

(The  traction  problem  has  also  been  considered  [21]).  We  want  to  show  that, 

A 

if  F,  g,  <t>  and  the  initial  history  y(0)  are  "nice",  then  the  initial  value 
problem  is  uniquely  solvable. 

D  k 

For  the  analysis,  we  introduce  the  following  function  spaces:  Wp'  ($1) 

denotes  the  usual  Sobolev  spaces,  and  we  write  wP'^ ,  Wp,,t  to  indicate  spaces 

of  vector-valued  or  symmetric  tensor-valued  func-tions.  p  is  assumed  to  lie 

between  3  and  6.  For  the  history  dependence  we  use  c^*m»  the  space  of 

bounded  continuous  functions  which  have  a  limit  at  -«•.  We  write  C^*m,  C?;im 

— b  =b 

for  vector-  and  tensor-valued  functions  and  C^im(X)  for  functions  taking 
values  in  the  Banach  space  X. 

The  following  describes  the  main  ideas  used  in  dealing  with  (3.1)- 
(3.5).  For  details,  the  reader  is  referred  to  [21]. 

1)  Of  course  the  basic  idea  is  to  treat  (3.4)  as  a  perturbation  of  the 
Navier-Stokes  equation.  The  theory  of  the  Navier-Stokes  equation  has 
been  developed  in  the  Euler ian  frame.  In  order  to  carry  over  Eulerian 
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results  to  Lagrangian  coordinates,  we  need  sufficient  smoothness  of  the 
transformation.  This  transformation,  however,  is  itself  one  of  the 
unknowns,  so  its  smooth ness  has  to  be  inferred  from  the  equation 
itself.  In  order  to  satisfy  this  consistency  requirement,  it  is 
advantageous  to  differentiate  (3.4)  with  respect  to  time.  We  also 

differentiate  (3.1)  and  (3.5)  twice  with  respect  to  time,  This  leads 

•  »•  _  • 

to  a  system  of  equations  for  y,  u  *  y,  a  *  y,  p  and  q  “  p. 

2)  Equation  (3.4)  and  the  first  time  derivative  of  (3.1)  can  be  used  to 

A 

express  u  and  p  in  terms  of  a  and  y*  This  involves  solving  an 
elliptic  system  129].  There  is  a  gain  of  regularity:  u  has  two  more 
derivatives  than  a . 

3)  The  boundary  and  incompressibility  constraints  for  a  can  be  reduced 
to  a  homogeneous  form  by  determining  an  appropriate  reference  function 
and  subtracting  it  from  a.  The  determination  of  this  reference 
function  again  involves  solving  an  elliptic  system. 

4)  A  projection  operator  is  used  to  eliminate  q.  This  leaves  an 
evolution  problem  for  y,  b  (the  reference  function  from  step  3)  and 
d  (a  new  variable  which  replaces  a). 

5)  In  order  to  deal  with  the  history  dependence,  we  now  regard  the  system 
from  step  4  as  an  evolution  problem  on  a  history  space.  In  doing  this, 
we  follow  the  following  recipe.  Suppose  you  have  an  equation  of  the 
form 

z  *  F(z,t)  ,  (3.6) 

A  <fe 

where  z(t)(s)  *  z(t+s),  s  e  (“•,<>].  'Bien  we  define  an  operator  F 

A  AAA 

which  maps  *  to  the  history  of  F:  F(z,t)(s)  *  F(Tsz,t+s),  where 
Tgz(t)(r)  -  z(t)(s+r)  -  z(t+s+r).  If  the  initial  history  satisfies  the 
equation,  (3.6)  can  be  written  in  the  form 
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I  A  A 

z  =  F(  z  1 1 ) 


(3.7) 


We  can  always  make  the  initial  history  satisfy  the  equation  by  adding 
an  appropriate  term  to  the  body  force. 

AAA 

6)  This  finally  leads  to  an  evolution  problem  for  y»  b,  d.  When  posed  on 
the  space  C^im(W^'4(fl)  *  W^'2(n)  x  £p(Q))  where  denotes  the 
subspace  of  divergence-free  vector  fields  with  zero  normal  component  on 
the  boundary,  this  problem  satisfies  the  assumptions  of  a  theorem  due 
to  Sobolevskii  [28]  on  abstract  quasilinear  parabolic  equations.  The 
essential  point  in  proving  this  is  of  course  the  fact  that  the  Stokes 
operator  generates  an  analytic  semigroup  in  iP  ([8],  [30]). 


4.  An  example  of  a  hyperbolic  model 

The  K-BKZ  model  [2],  [14]  is  motivated  by  an  analogy  with  finite 
elasticity.  The  constitutive  law  for  an  incompressible  elastic  material 
has  the  form 


„  -  +  ?W_  Y-1 

* "  PY  +  ar;  Yo 


23L  Y-\  y"1 

31  T  T0T 


(4.1) 


where  YQ  is  a  constant  tensor  and  W  is  a  scalar  function  of  1^  = 
trtYYp1)  and  X 2  =  tr(Y  *Yg).  Kaye  [14]  and  Bernstein,  Kearsley  and 
Zapas  substituted  the  following  for  a  viscoelastic  material 


-PY_1  +  ft.  a(t-T)[||-  y_1(t)  -  Y_1(t) 


(4.2) 


Y(T)Y-1(t)]dT 

with  I1  =  tr(Y(t)Y  Vt))  and  *  tr(Y  ^tJyfT)).  The  model  thus  assumes 
that  every  previous  state  of  the  material  is  like  a  temporary  equilibrium 
state  to  which  the  material  likes  to  revert;  the  influences  of  all  the 
previous  are  assumed  to  superpose  in  an  additive  fashion. 
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t 


We  assume  that  the  kernel  a  is  positive  and  smooth,  including 

t-T  =  0.  Let  F  denote  the  relative  deformation  gradient,  F*  = 

—  -  and  let  denote  the  entries  of  F-^:  F?  =  — - .  The 

3yp<T)  1  1  3/(1) 

equation  of  motion  can  be  written  in  the  form 


••i  3p  3?S  1  rt  . 

Py  *  “  “  -V  +  =■  J  _  a(t-T) 

acs  ay1  2  ^ 


.  32W  f  aV  3;q  3cs _  +  3y£ 

3f*3f^  3cq3CS  3yr(T)  3yp ( r )  3cq 


J_  (j£»)  dT  +  gi 

3;s  3yr (t )  3yp(r) 


(4. 


d.tt’L)  -  .  . 


3?' 


Hie  "hyperbolic"  character  of  (4.3)  is  guaranteed  by  a  strong  ellipticity 
condition,  which  has  the  same  form  as  in  elasticity 


r  32W  ,  — p— r i , i , j 

I - - -  +  KF ,  F  .  ]  X  XJy  y 

v  *  •>  i  j'  p  r 


3F13f;] 


P  r 


(4.. 


>  cm2  iwi2 


C  >  0 


for  large  enough  K.  This  condition  is  expressed  in  terms  of  F  and  has  a 
rather  indirect  form  in  terms  of  1^  and  I2.  However,  it  is  possible  to 
give  the  following  sufficient  condition:  (4.4)  holds  if  W  is  monotone  in 
both  1^  and  I2,  strictly  monotone  in  at  least  one  of  them,  and  W  is  a 
convex  function  of  /l1  and  /l2* 

In  [23],  I  proved  a  local  existence  theorem  for  (4.3)  posed  in  all  of 

2 

space,  assuming  that  (4.4)  holds.  The  analysis  proceeds  in  L  -type  spaces, 
i.e.  we  deal  with  solutions  for  which  £  -  £  ♦  0  at  infinity. 


I 


Equation  (4.3)  can  formally  be  written  as  a  non-delay  evolution  problem 
on  a  history  space.  In  this  case,  we  have  chosen  a  way  of  doing  this  which  is 
different  from  the  one  adopted  in  chapter  3.  For  each  t  >  0,  we  put 

A  A 

jj(t)(o)  m  y(ot).  That  is,  2.  only  contains  the  history  from  the  initial 
time  t  =  0  up  to  the  present  time,  not  the  history  for  t  <  0  (which  is  of 
course  considered  known). 

Again,  it  is  advantageous  to  differentiate  the  equation  with  respect  to 
time,  in  this  case  we  do  this  twice.  Lower  order  time  derivatives  can  be 

expressed  in  terms  of  higher  order  time  derivatives  by  solving  systems  which 

are  elliptic  in  the  sense  of  Agmon,  Douglis  and  Nirenberg  [1].  For  instance, 

(4.3)  can  be  regarded  as  a  nonlinear  elliptic  equation  for  y  and  p,  if  you 
presume  y  is  known.  One  then  looks  at  the  second  time  derivative  of 

(4.3) .  The  variable  p  can  be  eliminated  from  this  equation  using  the  Hodge 

projection.  One  finally  ends  up  with  an  evolution  problem  for  the  two 

variables  2  “  —  ”  £.)  and  *2"  ~  p  where  X  is  an  appropriately  chosen 

constant . 

In  this  system  the  leading  operator  is  elliptic  and  generates  a  quasi¬ 
contraction  semigroup.  An  existence  theory  for  systems  of  this  kind  is 
provided  by  Hughes,  Kato  and  Marsden  (11].  They  study  quasilinear  evolution 
equations  in  a  reflexive  Banach  space  X  of  the  form 

u  =  A(t,u)u  +  f(t,u)  ,  (4.5) 

where  A(t,u)  is  a  possibly  unbounded  linear  operator  and  f  is  a  bounded 
nonlinear  term.  The  essential  assumptions  of  the  theory  are  that  A 
generates  a  quasi-contraction  semigroup  (in  a  norm  which  is  allowed  to  depend 
on  t  and  u) ,  and  that  there  is  an  "elliptic"  operator  S(t,u)  such  that 
SAS-1  -  A  is  bounded.  Here  "elliptic"  means  that  S  is  a  bijection  from  an 
embedded  space  Y,  which  is  contained  in  the  domain  of  A  and  independent 
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of  t  and  u,  onto  X.  In  the  present  case,  A  turns  out  to  be  elliptic* 
and  we  can  simply  take  S  ■  A.  The  theorem  of  Hughes,  Kato  and  Marsden 
guarantees  the  local  existence  of  solutions  to  the  initial  value  problem,  if 
we  assume  sufficient  smoothness  of  the  data.  These  solutions  can  be  obtained 
by  the  iteration 

•n+1  _  n,  n+1  ,  n.  ,, 

u  *  A(t,u  )u  +  f(t,u  )  .  (4.6) 

The  most  interesting  question  from  the  rheologist's  point  of  view  is 
probably  whether  (4.4)  is  valid.  While  it  is  clear  that  (4.4)  would  hold  for 
small  deformations  in  any  reasonable  model,  the  global  validity  is  not  so 
clear.  Some  popular  rheological  models,  e.g.  W  =  I 1  or  I2,  always  satisfy 
(4.4).  For  other  models,  however,  (4.4)  fails  at  large  deformations.  This  is 
not  necessarily  bad.  If  (4.4)  fails,  the  evolutionary  character  of  the 
equations  can  be  lost  (cf.  also  [27]),  and  one  expects  something  strange  to 
happen  to  the  material.  In  fact,  strange  things  do  happen  at  high  shear 
rates,  and  loss  of  hyperbolicity  may  be  a  possible  explanation  for  these 
phenomena,  generally  known  as  "melt  fracture". 
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